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A  mathematical  model  for  the  performance  of  the  planar  solid  oxide  fuel  cell  (SOFC)  stack  is  described. 
The  model  considered  the  electric  contact  resistance  between  the  electrode  and  interconnect  rib,  the 
gas  transport  in  the  electrodes,  electronic  and  ionic  conductions  in  the  membrane-electrode  assembly 
and  the  electrochemical  reactions  at  the  gas-electrode-electrolyte  three  phase  boundaries.  The  model  is 
capable  of  describing  in  detail  the  rib  effect  on  the  gas  transport  and  the  current  distribution  in  the  fuel 
cell.  The  contact  resistance  is  found  to  be  an  important  factor  in  limiting  the  SOFC  performance.  Based 
on  the  interplay  of  the  concentration  and  ohmic  polarizations,  numerical  results  are  provided  for  the 
optimal  rib  widths  for  different  pitch  sizes  and  different  area  specific  contact  resistance  (ASRCOntact)-  The 
optimal  rib  width  is  found  to  be  linear  to  the  pitch  width  for  a  given  ASRCOntact  and  the  parameters  for  the 
linearity  are  given.  The  parameters  are  little  affected  by  the  hydrogen  concentration  and  the  thickness, 
porosity  or  conductivity  of  the  cathode.  The  influence  of  the  cathode  thickness  on  the  SOFC  performance 
is  also  examined.  Contrary  to  the  common  belief  on  the  thin  cathode  (~50  |mm),  thicker  cathode  layer 
(100-300  |jim)  is  beneficial  to  the  SOFC  stack  performance. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cell  (SOFC)  offers  great  potential  for  high 
efficiency,  environmentally  benign  electric  energy  generation  by 
converting  chemical  energy  to  electricity  directly  without  pass¬ 
ing  through  the  combustion  process  and  Carnot  circle.  There  is 
a  great  deal  of  research  on  solid  oxide  fuel  cells  generally,  espe¬ 
cially  on  new  material  systems  and  on  the  optimizations  of  the 
microstructures  to  improve  the  electrochemical  performance  of  the 
cell,  or  membrane-electrode  assembly  (MEA).  Planar- type  designs 
have  received  much  attention  because  they  potentially  offer  higher 
power  density  relative  to  the  tubular-type  SOFC,  which  is  ascribed 
to  the  low  electrical  resistance  due  to  shorter  current  paths.  In 
particular,  a  very  thin  electrolyte  film  can  be  used  in  an  anode- 
supported  SOFC,  drastically  reducing  the  ohmic  resistance  and 
enabling  operation  at  intermediate  temperature  [1-3].  However, 
relatively  little  of  the  literature  concentrates  on  such  critical  tech¬ 
nical  issues  as  the  system-level  layout  of  the  layers  and  stacks  [4-6]. 

The  most  popular  and  reliable  MEA  materials  for  SOFCs  com¬ 
prise  Y203 -stabilized  zirconia  (YSZ)  as  the  electrolyte  and  a 
constituent  in  the  composite  anode  and  cathode,  Ni  for  the  anode, 
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and  Sr-doped  LaMn03  (LSM)  for  the  cathode.  Under  ideal  condi¬ 
tions  with  pure  hydrogen  fuel  and  low  fuel  utilization,  numerous 
groups  have  reported  over  the  years  that  single  cells  with  optimized 
materials  can  yield  power  densities  on  the  order  of  1 W  cm-2  or 
higher  at  the  operating  temperature  of  700-800  °C  [2,3,7-12].  For 
practical  applications,  however,  multiple  cells  must  be  connected 
in  series  by  interconnect  material  to  form  a  fuel  cell  stack.  The 
recently  reported  best  performance  of  SOFC  stacks  with  metallic 
bipolar  plate  current  collectors  [13,14]  is  substantially  below  that 
for  the  single  cell  [2,3].  The  lowered  cell  performance  in  the  stack 
may  be  attributed  to  the  reduced  gas  transport  and  the  increased 
ohmic  polarization  in  the  stack  setting.  When  assembling  the  SOFC 
stack,  the  connecting  materials  (bipolar  plates)  are  required.  Small 
channels,  which  are  often  formed  in  the  interconnect  material,  are 
commonly  used  to  carry  the  fuel  and  air  gas  flow.  The  “ribs”,  which 
separate  and  define  the  flow  channels,  make  direct  contact  with  the 
electrodes.  For  technical  reason,  the  channels  and  the  ribs  cannot  be 
too  small  [15]  and  the  contact  resistance  between  the  ribs  and  the 
electrodes  is  part  of  the  intrinsic  property  of  the  interconnect  mate¬ 
rial  [16].  In  designing  the  layer  architecture,  there  is  a  tradeoff  that 
must  be  considered  between  the  rib  and  channel  sizes.  On  one  hand, 
wider  ribs  and  ribs  covering  bigger  fraction  of  the  cell  area  may 
reduce  the  interface  resistance  to  current  flow  by  increasing  the 
electrode-interconnect  contact  area  and  reducing  the  current  path 
through  the  possibly  high  resistance  electrode  material.  Hence, 
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such  ribs  will  give  a  better  conduction  of  the  electrical  current  and 
reduce  ohmic  losses.  On  the  other  hand,  the  chemical  species  do 
not  diffuse  as  well  underneath  wide  ribs.  Narrow  ribs  are  needed  to 
facilitate  more  uniform  distribution  of  reactive  gases  across  the  area 
of  the  electrolyte  surface  and  thus  to  promote  electrochemical  per¬ 
formance.  The  implication  of  the  tradeoff  to  the  cell  performance 
can  be  very  significant  and  a  designer  must  decide  how  to  configure 
a  channel  network  to  optimize  the  cell  performance. 

There  have  been  a  few  pioneering  works  on  the  rib  design  opti¬ 
mization.  Ferguson  et  al.  [17]  gave  a  numerical  example  for  the  rib 
size  effect  on  the  performance  of  a  specific  electrolyte-supported 
SOFC.  Lin  et  al.  [18]  provided  a  phenomenological  model  and  ana¬ 
lytical  expressions  to  estimate  the  rib  effects  on  the  concentration 
and  ohmic  polarizations  of  anode-supported  SOFC  stacks.  Ji  et  al. 
[19]  showed  that  the  terminal  output  of  a  SOFC  stack  depended 
strongly  on  the  contact  resistance.  Jeon  et  al.  [15]  proposed  a 
detailed  microstructural  model  and  examined  systematically  the 
influence  of  the  rib  width,  pitch  width  and  area  specific  contact 
resistance  (ASRCOntact)  on  the  stack-cell  performance.  The  work  by 
Jeon  et  al.  provided  very  useful  information  in  revealing  the  details 
of  the  complex  interaction  of  competing  factors  in  SOFC  operations. 
Unfortunately,  their  optimization  results  appear  to  be  dependent 
on  a  large  number  of  design  and  operational  parameters  and  are 
difficult  to  generalize  and  to  provide  easy-to-use  design  guidance 
to  the  broad  fuel  cell  engineering  society.  Moreover,  their  treatment 
of  the  effective  in-plane  area  specific  resistance  appears  to  be  over¬ 
simplified  and  may  affect  the  optimization  result  for  the  cathode 
current  collect  layer  thickness. 

In  the  present  study,  the  cell  performance  of  an  SOFC  stack  is 
analyzed  with  a  mathematical  model  capable  of  describing  the 
details  of  gas  transport  in  the  electrodes,  electronic  and  ionic  con¬ 
ductions  in  the  MEA,  the  electrochemical  reactions  at  the  three 
phase  boundaries  (TPBs)  and  the  current  distributions  in  the  stack¬ 
cell.  The  optimal  rib  designs  for  different  combinations  of  pitch 
widths  and  ASRCOntact  are  obtained  and  parameterized  to  provide 
an  easy-to-use  guidance  for  optimizing  the  rib-channel  layout.  The 
applicability  of  the  results  to  the  stack  operations  is  validated  by 
sensitivity  analysis  for  a  range  of  hydrogen  concentrations.  The  gen¬ 
erality  of  the  results  for  different  fuel  cells  is  illustrated  by  testing 
on  electrodes  with  different  thickness,  different  porosity  and  dif¬ 
ferent  conductivities.  The  influence  of  the  cathode  thickness  on  the 
performance  of  a  fuel  cell  stack  is  also  examined. 


2.  Method 

We  limit  ourselves  to  the  discussion  of  anode-supported  pla¬ 
nar  SOFC  stacks  with  co-flow  or  counter  flow  designs  to  simplify 
the  two-dimensional  (2D)  model  setup.  The  modeling  results  for 
the  rib  size  optimization  might  be  similar  for  the  cross-flow  design 
as  the  obtained  results  to  be  discussed  below  are  only  sensitive  to 
the  rib-electrode  contact  resistance.  A  schematic  cross-section  of  a 
repeating  unit  is  sketched  in  Fig.  1.  The  fuel  and  oxidant  gases  are 
supplied  to  the  cell  through  the  channel/electrode  boundary  and 
the  current  is  collected  through  the  rib.  As  the  TPBs  are  effective 
only  at  the  region  very  close  to  the  electrode/electrolyte  interfaces, 
the  TPBs  are  assumed  to  be  located  only  at  the  electrode/electrolyte 
interfaces  and  there  are  no  electrochemical  reactions  inside  the 
electrodes  in  our  numerical  model.  Though  metallic  interconnects 
with  high  conductivities  are  the  desired  and  frequently  used  mate¬ 
rials  for  the  bipolar  current  collectors,  the  contact  resistance  at  the 
rib  and  electrode  boundary  may  be  substantial  due  to  the  forma¬ 
tion  of  oxide  scales.  Many  efforts  are  devoted  to  the  reduction  of  the 
area  specific  contact  resistance  of  low  cost  alloy  interconnects,  such 
as  Ni-,  Fe-,  and  Co-based  superalloys,  Cr-based  alloys,  and  stain¬ 
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Fig.  1.  Schematic  cross-section  of  a  repeating  unit  in  a  planar  SOFC  stack. 


less  steel,  to  the  desired  level  of  25-50  m^  cm2  [16,20,21].  At  the 
desired  level  of  contact  resistance,  the  fuel  cell  performance  may 
still  be  significantly  affected  and  the  design  optimization  is  impor¬ 
tant.  The  general  approach  in  discussing  the  rib  effects  is  to  map  out 
the  transportation  process  in  the  case  of  rib  presence,  starting  from 
the  known  parameters  when  ribs  are  absent.  The  gas  transportation 
path  determines  the  concentration  polarization  and  electrochem¬ 
ical  performance,  and  the  electrical  path  determines  the  ohmic 
polarization.  The  overall  cell  performance  as  well  as  the  optimal 
rib  design  can  then  be  determined. 

2.1.  Physical  model 

Due  to  symmetry  reason,  only  half  of  the  repeating  unit  shown 
in  Fig.  1  needs  to  be  examined.  Furthermore,  the  gas  channels  and 
the  bipolar  plate  can  be  replaced  with  the  proper  boundary  condi¬ 
tions  for  gas  transport  and  electric  potential.  The  contact  resistance 
is  set  on  the  boundary  of  interconnect  ribs  and  electrodes.  With 
good  connection  between  interconnect  and  electrode  materials,  the 
contact  resistance  is  dominated  by  the  oxide  scale  of  the  intercon¬ 
nect  [16,20,21  ]  and  is  assumed  here  to  be  uniform  in  the  rib  region. 
Moreover,  the  ASRCOntact  on  the  anode  side  and  on  the  cathode  side 
are  set  to  be  the  same  in  our  model  and  is  equivalent  to  the  aver¬ 
age  ASRcontact  of  the  anode  and  cathode  sides  as  the  total  currents 
for  the  anode  and  cathode  ribs  are  identical.  In  our  calculations  for 
the  thin  cathode  layer  stacks,  it  is  noticed  that  part  of  the  cath¬ 
ode  area  underneath  the  interconnect  rib  may  run  out  of  gas  at 
high  current  density.  The  gas  depletion  zone  may  cause  numerical 
instability  in  simulating  the  gas  transport  process  as  the  bound¬ 
ary  condition  derived  from  the  electrochemical  reaction  may  cause 
negative  oxygen  concentration  there.  Similar  numerical  instability 
has  been  encountered  by  other  researchers  [15].  To  overcome  the 
numerical  difficulty,  the  gas  depletion  area  may  be  identified  with  a 
cutoff  gas  pressure  (0.5%pO2 )  and  treated  as  a  pure  electric  conduct¬ 
ing  area.  The  cutoff  pressure  has  been  verified  to  cause  negligible 
error  in  the  result.  For  example,  the  average  output  current  density 
for  the  baseline  case  with  rib  width  of  0.8  mm  to  be  described  in 
Section  3.1  is  0.5143  and  0.5145  A  cm-2  for  the  cutoff  gas  pressure 
of  0.5%  and  0.1%,  respectively.  Except  the  pure  conducting  area,  the 
cathode  is  treated  as  the  usual  gas  diffusion  and  electric  conduc¬ 
tion  area.  As  a  result,  the  overall  physical  model  is  reduced  to  that 
shown  in  Fig.  2.  For  easy  description  of  boundary  conditions  in  later 
part  of  the  paper,  each  boundary  line  segment  in  our  2D  model  is 
assigned  with  an  alphabetic  name  as  indicated  in  Fig.  2. 
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Fig.  2.  Sketch  of  the  physical  model  and  subdomains  and  boundary  designations. 


The  overall  cell  performance  can  be  described  by  the  operating 
cell  potential  and  the  output  current.  The  operating  cell  potential 
(Vcell)  can  be  formally  expressed  as: 


Vceii  —  E0  PASR;a  Pconc;a  Pact;a  Pohm;a  Pohm;el  Pohm;c 

—  Pact;c  —  Pconc;c  ~  PASR;c  (1) 

where  E0  is  the  Nernst  potential,  r]conc-a  and  r]conc-c  are  respectively 
the  concentration  overpotentials  at  the  anode  and  the  cathode  due 
to  the  gas  species  diffusion  resistance,  rjac t;a  and  qac t;c  the  corre¬ 
sponding  activation  overpotentials,  i)ohm.at  q0 hm;ei  and  p0hm;c  the 
ohmic  overpotentials  in  the  anode,  electrolyte  and  cathode,  rjASR-,a 
and  ?7asr;c  the  anode  and  cathode  interface  overpotentials  due  to 
the  contact  resistance  at  material  boundaries.  All  these  overpo¬ 
tentials  are  functions  of  the  current  density,  j.  Models  for  each  of 
the  overpotentials  will  be  described  in  subsequent  sections  in  the 
paper.  As  Vcell  is  an  operating  parameter,  Vcell  is  used  as  an  appropri¬ 
ate  boundary  condition  for  solving  the  current  density  distribution 
in  the  physical  model. 

The  average  output  current  density  of  the  SOFC  stack  can  be 
calculated  by 


J  = 


^  /‘“pitch 

tfpitch  Jx= o 


jydX 


(2) 


where  dpitch  is  the  pitch  width  (the  sum  of  the  rib  width,  drib,  and 
the  channel  width,  dchannei).  The  x-  and  y-axis  are  along  the  hori¬ 
zontal  and  vertical  directions  in  the  2D  model,  respectively.  jy  is  the 
y  component  of  the  current  density  flux  vector. 


2.2.  Gas  transport  modeling 


The  gas  flow  transports  in  the  anode  and  cathode  are  governed 
by  the  momentum  and  mass  balance  equations. 


2.2.1.  Momentum  balance  equation 

In  the  porous  structure  of  electrodes,  the  gas  flow  is  mainly 
driven  by  the  pressure  gradient.  The  Darcy’s  law  and  the  continuity 
equation  were  used  to  simulate  the  flow  in  electrodes 


u  =  -  - Vp  (3a) 

p 

V  •  pu  =  0  (3b) 


where  k  denotes  the  permeability  of  the  porous  media  (m2),  q 
the  fluid  viscosity  (Pas),  p  the  pressure  (Pa),  u  the  velocity  vec¬ 
tor  (ms-1)  and  p  is  the  density  of  the  fluid  (kg  m-3 ).  The  density  of 
fuel  (pf)  and  air  (pair)  are  calculated  by  the  gases  molar  fractions 


Pf  = 


P(xh2Mh2  +  *h2oMh2o) 
RT 


Pair  = 


p(xo2M02  +xN2MN2) 

RT 


(4a) 

(4b) 


where  M,-  is  the  molar  mass  of  species  i,  Xj  the  molar  fraction  of 
species  j,  R  the  universal  gas  constant  and  T  is  the  temperature. 


The  anode  concentration  overpotential  is  related 
to  the  partial  pressures  of  fuel  species  by  pConc;a  = 
(i?r/2F)ln((pH2,f/PH2o,f)(PH2o,TPB/PH2,TPB)),  where  F is  the  Faraday 
constant.  The  anode  concentration  overpotential  depends  on  the 
x  coordinate  of  the  TPB  in  anode.  Similarly,  the  cathode  concen¬ 
tration  overpotential  is  related  to  the  oxygen  partial  pressures  by 
Pconc;c  =  (P’P/4F)ln(po2?air/po2jpB)  and  depends  on  the  x  coor¬ 
dinate  of  the  TPB  in  cathode.  In  this  article,  subscript  a,  c,  el,  I, 
air,  f  and  TPB  refers  to  anode,  cathode,  electrolyte,  interconnect, 
air,  fuel  and  TPB,  respectively.  As  the  Nernst  potential  is  given 
as 


where  Eg  is  the  Nernst  potential  when  the  partial  pressure  of 
H2,  H20  and  02  are  all  at  1  atm  (and  Eg  =  1.01  V  at  E=973  K  or 
700  °C),  the  last  two  terms  of  Eq.  (5)  can  be  combined  with  qCOnc;a 
and  ?7conc;ct  respectively,  to  yield  the  anode  and  cathode  balance 
potentials: 


f  P n2,a  \ 
VPH20,a  ) 


TPB 


0  RTi 

fi  =  -4 f'n 


f  Po2,cItpb\ 
y  1  atm  j 


(6a) 

(6b) 


Eq.  (1 )  may  then  be  rewritten  as 

Vceii  =  Eq  —  ?7ASR;a  —  Pa  —  Pact;a  —  Pohm;a  —  Pohm;el  —  Pohm;c  —  Pact;c 

-^c-'Msr*  O') 


2.2.2.  Mass  balance  equation  (Maxwell-Stefan  equation) 

In  the  SOFC  porous  electrodes,  the  relative  concentrations  of  gas 
species  are  of  the  same  order  of  magnitude.  The  Maxwell-Stefan 
equation  is  suitable  for  calculating  the  gas  transport  processes.  The 
governing  equation  is 


V  •  ^pa)iU  -  pOJi'y'pij(VXj  +  (Xj  - 


=  0 


!>>•= 


(7a) 


(7b) 


where  cO[  is  the  mass  fraction  of  species  i. 

D  is  the  multi-component  Fick  diffusivity.  It  is  related  to  the 
multi-component  Maxwell-Stefan  diffusivity  (D)  by  the  expres¬ 
sions  [22] 


coi 

Dl1  =  ^Dl2 
D12=D21=-^D12 


*1*2 


°22  =  ^ ~2Du 

The  multi-component  Maxwell-Stefan  diffusivity  can  be  calcu¬ 
lated  by  the  expression  [23] 

i1/2 


(8a) 

(8b) 

(8c) 


Dii  =  3.198  x  10 


i-8 


7-1.75 


p(v,1/3+v;/3) 


i  i 

Wi  +  Wi 


(9) 
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where  vz  is  the  diffusion  volume  for  species  i.  In  porous  media,  the 
effective  diffusivity  is  reduced  to  [24] 

Df  =  eTDij  (10) 

where  £  is  the  porosity  and  r  the  tortuosity  of  the  porous  media. 

2.3.  Electric  current  transport 

The  electrical  current  is  determined  by  the  charge  balance 


the  partial  pressure  and  the  activation  polarization  (overpotential) 
by  the  Butler-Volmer  equation.  However,  different  parameteriza¬ 
tion  forms  have  been  proposed  based  on  different  experimental 
measurements.  As  long  as  an  expression  is  able  to  represent  the 
experimental  data  reasonably  well,  the  specific  form  is  not  critical 
due  to  the  use  of  fitting  parameters. 

According  to  the  empirical  forms  suggested  before  [15,25-27], 
the  Butler-Volmer  equation  used  here  for  the  anode  side  TPB  is 
expressed  as 


forpH2o,TPB  <  14,  000  Pa 
forpH2o,TPB  >  14,  000  Pa 


(12a) 


equation.  Electronic  charge  transfer  equation  was  used  to  obtain 
the  current  distribution  in  the  anode  and  cathode  area 

-V-(oeWe)  =  0  (11a) 

where  cre  is  the  electronic  conductivity  of  the  electrode  material 
and  Ve  the  electric  potential  in  the  electrode  and  -crevVe  the  flux 
vector  of  the  electronic  current  density.  The  difference  between 
the  electric  potentials  along  the  electronic  current  flux  path  at  the 
electrode-rib  boundary  and  the  TPB-electrode  boundary  is  the  elec¬ 
trode  ohmic  polarization  (p0hm;a  or  ^ohm;c)  and  is  generally  current 
path  dependent. 

Ionic  charge  transfer  equation  was  used  in  the  electrolyte 


-V.(aiWi)  =  0 


(Hb) 


where  ax  is  the  ionic  conductivity  of  the  electrolyte,  Vx  the  elec¬ 
tric  potential  in  the  electrolyte  and  -crxvVx  is  the  flux  vector  of  the 
ionic  current  density.  The  difference  between  the  electric  potentials 
along  the  ionic  current  flux  path  at  the  electrolyte-anode  bound¬ 
ary  and  the  electrolyte-cathode  boundary  is  the  electrolyte  ohmic 
polarization  (p0hm;ei)  and  is  also  current  path  dependent. 

The  electric  potential  loss  inside  interconnect  plate  is  assumed 
to  be  negligible  due  to  high  conductivity  of  the  metallic  material. 
The  local  current  densities  cross  the  interconnect/anode  (j\^a)  and 
the  cathode/interconnect  (jc^i)  interfaces  are  determined  by  the 
associated  electric  potential  changes,  or  the  interface  overpoten¬ 
tials: 


J I-»a  = 


Jc— >1  = 


^e,I/a  ^e,a/I 
ASRcontact 
Ve,c/I  ~  Ve,I/c 
ASRcontact 


^ASR;a 

ASRcontact 

(lie) 

*7aSR;c 

ASRcontact 

(lid) 

where  VeI/a  and  Ve  are  respectively  the  interconnect  and  anode 
electric  potentials  at  the  anode-interconnect  boundary,  Ve  >c/i  and 
VeI/c  the  cathode  and  interconnect  electric  potentials  at  the 
cathode-interconnect  boundary. 


where  ja  is  the  current  density  and  j0  a  the  exchange  current  density 
for  the  anode  TPB.  The  cutoff  pressure  of  14,000  Pa  for  Ph2o,tpb  is 
used  to  account  for  the  catalytic  effect  of  water  vapor  [28]. 

The  Butler-Volmer  equation  used  for  the  cathode  side  TPB  is 
[26,29,30] 


where  jo, c  is  the  exchange  current  density  for  the  cathode  TPB. 

The  activation  polarizations,  rjact,a  and  rjact,c,  are  related  to  the 
electric  and  balance  potentials  by 


*7act,a  —  ^e,a/el  -  ^i,el/a  -  ha  (13a) 

hact,c  —  ^i,el/c  —  ^e,c/el  —  he  (13b) 

where  Ve  a/el  and  Viel/a  are  respectively  the  anode  and  electrolyte 
electric  potentials  at  the  anode-electrolyte  boundary,  Vx  el/c  and 
Ve  c/el  the  electrolyte  and  cathode  electric  potentials  at  the  cathode- 
electrolyte  boundary. 


2.5.  Boundary  conditions  (BCs) 

As  described  above,  the  governing  partial  differential  equations 
(PDEs)  are  Darcy’s  law  and  Maxwell-Stefan  equation  for  the  gas 
diffusions  in  the  anode  and  cathode,  Ohm’s  law  for  electrical  cur¬ 
rent  conduction  in  the  anode,  cathode  and  electrolyte.  The  settings 
for  the  boundary  conditions  required  for  solving  these  PDEs  are 
summarized  in  Tables  1  and  2. 

The  Butler-Volmer  equations  are  coupled  with  the  governing 
equations  and  provide  the  current  densities  at  the  TPB  boundaries 
for  given  gas  concentrations  and  activation  overpotentials.  The  cur¬ 
rent  densities  then  serve  as  the  boundary  conditions  to  solve  the 
governing  equations,  resulting  in  new  TPB  gas  concentrations  and 
activation  overpotentials.  The  iterative  process  continues  until  a 
convergent  solution  is  found. 


2.4.  The  current  density  at  the  TPB 


2.6.  Basic  model  parameters 


The  current  density  at  the  TPB  is  proportional  to  the  electro¬ 
chemical  reaction  rates.  It  is  related  to  the  exchange  current  density, 


Unless  stated  otherwise,  the  following  model  parameters  are 
used.  Cell  operates  at  T=973I<  (700  °C)  with  both  fuel  and  air  at 


Table  1 

Boundary  settings  for  Darcy’s  law  and  the  Maxwell-Stefan  equation 


Boundary 

a 

h 

0 

n 

All  others 

Darcy’s 

law 

Maxwell-Stefan 

equation 

BC  type 

BC 

BC  type 

BC 

Pressure  condition 

1  atm 

(02)  mass  fraction 

0.23  (air) 

(H2)  mass  fraction 

< 

Velocity  inflow/outflow 
J^(Mh2o-MH2) 

(02)  flux 

-jcM02 

4  F 

Velocity  inflow/outflow 

4ffeM°2 

(H2)  flux 

~jaMH2 

2  F 

Insulation/symmetry 

Insulation/symmetry 
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Table  2 

Boundary  settings  for  the  electronic  and  ionic  charge  transfer  equations 


Boundary 

b,  c 

g 

0 

n 

m 

All  others 

Electronic 

Ionic 

BC  type 

BC 

BC 

Reference  potential 
Vcen-Eg 

Reference  potential 

0 

Inward  current  flow 
jc 
jc 

Inward  current  flow 

-Ja 

ja 

Inward  current  flow 

0 

0 

Electric  insulation 

pressure  of  1  atm.  The  molar  fraction  of  hydrogen  in  the  fuel  is 
0.90.  The  output  cell  voltage  is  fixed  at  0.7  V,  a  balanced  choice  for 
cell  efficiency  and  power  density  [24].  The  viscosities  for  fuel  and 
air  are  respectively  2.8  x  10-5  and  4  x  10-5  Pas  [24,31].  The  diffu¬ 
sion  volumes  for  H2,  H20,  02  and  N2  are  7.07  x  10-6,  12.7  x  10-6, 
16.6  x  10-6  and  17.9  x  10-6  m3  mol-1 ,  respectively  [32].  Anode  per¬ 
meability  and  cathode  permeability  are  respectively  1.7  x  10-10 
and  1  x  10-13  m2  [31,33].  The  thicknesses  for  anode,  cathode  and 
electrolyte  layers  are  750,  50  and  10  [xm,  respectively.  Anode 
porosity  is  0.38  and  cathode  porosity  is  0.3.  Both  anode  and  cath¬ 
ode  tortuosities  are  3.  The  anode  and  cathode  exchange  current 
densities  are  4280  and  1070 Am-2,  respectively.  The  expres¬ 
sions  for  anode,  cathode  and  electrolyte  conductivities  (sm-1) 
are  aa  =  3.356  x  104  exp(1392/T),  orc  =  1.223  x  104  exp(-600/T)  and 
<rel  =  3.34  x  104exp(— 10,30 0/T),  respectively  [17,34].  The  fuel  cell 
described  with  these  parameters  is  sometime  referred  as  the  refer¬ 
ence  cell  for  convenience. 

2.7.  Numerical  solutions 

The  finite  element  commercial  software  COMSOL 
MULTIPHYSICS®  Version  3.4  [35]  was  used  in  the  present  study 
to  solve  the  PDEs  with  the  appropriate  boundary  settings.  The 
COMSOL  stationary  nonlinear  solver  uses  an  affine  invariant  form 
of  the  damped  Newton  method  [36]  to  solve  the  discretized  PDEs 


with  a  relative  tolerance  of  1  x  10-6.  Fine  triangle  meshes  were 
used  and  verified  to  provide  accurate  results  in  comparison  with 
the  very  fine  mesh  of  a  maximum  element  size  of  1  x  10-5  m 
for  boundaries  a,  b  and  c  and  0.5  x  10-6  m  for  region  around  the 
a  and  b  contact  area.  The  output  current  densities  for  various 
combinations  of  cell  designs  and  operating  parameters  were 
calculated.  As  many  calculations  were  involved,  the  COMSOL  Script 
1.2  [35]  was  employed  to  carry  out  the  algorithm  more  efficiently. 
The  examined  hydrogen  molar  concentrations  in  the  fuel  were 
0.97,  0.9,  0.6  and  0.29.  These  values  may  be  used  to  represent 
hydrogen  depletion  as  hydrogen  flow  down  the  channel.  The  rib 
effects  on  the  stack  operations  were  evaluated  with  the  specific 
area  contact  resistances  of  1  x  10-5,  0.01,  0.025,  0.05,  0.075,  and 
0.1  £2  cm2  for  the  interconnect-electrode  contacts.  Unless  stated 
otherwise,  the  specific  area  contact  resistance  of  0.05  £2  cm2  is 
used  for  the  interconnect-electrode  contacts.  The  pitch  width 
was  varied  from  0.5  to  5  mm  in  step  of  0.5  mm.  For  a  given  pitch 
width,  the  optimized  rib  width  was  found  by  varying  the  channel 
width  to  produce  the  maximum  output  power.  The  sensitivity  of 
the  optimization  results  on  various  factors  such  as  the  cathode 
thickness,  porosity  and  conductivity  were  also  tested.  The  effects 
of  cathode  thickness  on  the  cell  performance  were  examined 
by  varying  the  cathode  thickness  from  50  to  400  [xm.  Single  cell 
performance  used  below  are  calculated  using  the  model  without 
pure  conducting  area  and  no  contact  resistance,  and  all  the  pitch 
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Fig.  3.  Distributions  of  hydrogen  and  oxygen  molar  fractions  (*h2  ,  *o2  )• 
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Fig.  4.  The  distributions  of  the  anode  and  cathode  electric  potentials. 
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(channel  and  rib)  boundaries  were  set  as  gas  inlet/outlet  and 
electric  conductive  boundaries. 

3.  Results  and  discussion 

3.1.  Performances  of  cells  with  ribs 

We  begin  by  showing  specific  examples  of  the  effects  of  rib 
designs  on  the  fuel  cell  performances.  The  pitch  and  rib  widths  for 
the  stack  cell  are  2  and  0.8  mm,  respectively.  With  these  settings, 
the  mean  output  current  density  for  the  reference  cell  was  calcu¬ 
lated  to  be  0.51  A  cm-2,  a  significant  reduction  from  the  result  of 
1.4  A  cm-2  for  the  single  cell  without  the  ribs. 

Fig.  3  shows  the  hydrogen  and  oxygen  molar  fraction  distri¬ 
butions  over  the  anode  and  cathode,  respectively.  The  H2  molar 
fraction  varies  over  the  range  of 0.84-0.90.  The  maximum  H2  molar 
fraction  is  naturally  at  the  channel-anode  boundary,  and  the  min¬ 
imum  is  at  the  TPB  underneath  the  channel.  The  rib  effect  on  the 
anode  gas  transport  is  quite  limited.  While  for  cathode,  an  oxygen 
depletion  zone  with  a  width  of  0.46  mm  (23%  of  the  pitch  width) 
was  identified.  This  oxygen  depletion  zone  alone  would  contribute 
to  23%  reduction  of  the  cell  output.  As  shown  in  Fig.  3,  the  oxy¬ 
gen  concentration  is  very  low  in  the  diffusion  region  neighboring 
the  depletion  zone,  further  contributing  to  the  reduction  of  the 
stack  cell  performance.  The  diffusion  of  oxygen  is  much  easier  than 
hydrogen  in  the  vertical  direction  due  to  the  relative  short  diffusion 
path.  In  the  horizontal  direction,  however,  the  oxygen  diffusion  is 
limited  by  the  relative  narrower  alleyway  which  is  direct  propor¬ 
tion  to  the  cathode  thickness,  and  may  be  depleted  easily  when  the 
output  current  is  not  small. 

Fig.  4  shows  the  electric  potential  distribution  in  the  anode  and 
cathode.  As  stated  above  (see  Table  3),  the  reference  potential  at 
the  rib-anode  boundary  was  0  V  and  at  the  rib-cathode  bound- 
ary  Vce\\  -Eq  =  -0.31  V.  As  can  be  seen,  the  potential  loss  due  to 


Table  3 

Parameters  for  the  optimal  rib  width,  Eq.  (14),  for  various  ASRCOntact  (£2  cm2) 


the  anode-rib  contact  resistance  was  about  64  mV  and  the  anode 
side  potential  variation  is  very  small.  Notable  potential  variation  is 
observed  in  the  cathode.  The  potential  loss  due  to  the  cathode-rib 
contact  resistance  varied  from  60  to  72  mV.  The  cathode  potential 
varied  from  -0.250  to  -0.223  V  and  the  maximum  potential  change 
inside  cathode  was  27  mV.  The  ohmic  overpotentials  due  to  the  rib- 
electrode  contacts  are  significant  contributors  to  the  reduction  of 
the  stack  cell  performance.  The  cathode  ohmic  overpotential  is  also 
big  enough  to  influence  the  cell  output. 

As  the  cathode  conductivity  is  relatively  high,  large  cathode 
ohmic  overpotential  is  somewhat  unexpected.  It  becomes  clear 
when  the  distribution  of  the  current  flux  is  displayed.  Fig.  5  shows 
the  distributions  of  the  current  densities  and  flux  vectors  in  the 
anode  and  cathode  areas.  The  total  current  densities  in  the  cathode 
are  much  larger  than  the  stack  cell  output  current  density.  For  the 
current  generated  in  the  TPB  zone  to  reach  the  cathode  intercon¬ 
nect,  most  of  the  current  has  to  flow  along  the  horizontal  direction 
with  a  narrow  cross-section  defined  by  the  cathode  thickness  and  a 
long  path  length  determined  by  the  channel  width.  Consequently, 
the  cathode  ohmic  overpotential  for  the  stack  cell  (~10  mV)  is  sev¬ 
eral  orders  of  magnitude  higher  than  that  for  a  single  cell  with  Pt  or 
Ag  mesh  as  current  collector  (~0.04mV).  As  noted  by  passing,  the 
current  density  around  the  cathode,  channel  and  rib  boundary  zone 
is  particularly  high  and  may  be  a  significant  heat  source  and  cause 
mechanical  failures.  More  attentions  should  be  paid  to  the  area  in 
the  process  of  evaluating  the  mechanical  performance  of  the  fuel 
cell  stack. 
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Fig.  5.  The  current  density  distributions:  (a)  the  current  density  and  current  flux  vector  in  the  electrodes;  (b)  the  y-component  of  the  current  flux  vector  in  the  fuel  cell. 


220 


S.  Liu  et  al.  /  Journal  of  Power  Sources  183  (2008)  214-225 


x  (mm) 


c 

o 


o 

CD 


_cd 

o 


x  (mm) 


c 

o 

o 

CD 


CD 

O 


Fig.  6.  Distributions  of  the  TPB  current  density,  H2  and  02  molar  fractions:  (a) 
dpitch  2  mm  and  0,8  mm,  (b)  dpitch  5  mm  and  2  mm. 


Fig.  6  shows  the  hydrogen  molar  fraction  and  current  density 
distributions  at  the  anode  TPB  and  the  oxygen  molar  fraction  dis¬ 
tribution  at  the  cathode  TPB  for  two  stack  cell  designs  of:  (a) 
c/ pitch  =  2  mm,  dri5  =  0.8  mm  and  (b)  dpitch  =  5  mm,  drib  =  2mm.  For 
the  stack  cell  with  dpitch  =  2  mm  and  drib  =  0.8  mm,  the  current  den¬ 
sity  in  region  underneath  the  channel  (x<  1.2  mm)  is  quite  uniform 


and  varies  between  0.70  and  0.76  A  cm-2  with  a  maximum  at 
x  =  1.13  mm.  The  current  density  for  the  TPB  underneath  the  rib 
(x>1.2mm)  decreases  rapidly  with  increasing  x  and  is  negligible 
for  x  >  1.6  mm.  With  increasing  x  for  the  TPB  underneath  the  rib,  the 
H2  molar  fraction  increases  slightly,  while  the  O2  molar  fraction  of 
oxygen  decreases  rapidly.  Clearly,  the  TPB  current  density  distri¬ 
bution  is  mainly  determined  by  the  oxygen  concentration.  For  the 
stack  cell  design  of  dpitch  =  5  mm  and  drib  =  2  mm  that  has  the  same 
rib  to  pitch  width  ratio  as  the  former  design  of  dpitch  =  2  mm  and 
drib  =  0.8mm,  the  cell  output  current  is  reduced  to  0.39  A  cm-2,  a 
deduction  of  24%  from  the  cell  output  current  density  of  0.51  A  cm-2 
for  the  former  design.  The  performance  reduction  is  mainly  due 
to  the  increased  cathode  current  path  length  and  the  increased 
cathode  ohmic  overpotential  for  the  increased  channel  width,  as 
discussed  above.  With  the  wider  pitch  design,  the  TPB  current  den¬ 
sity  underneath  the  channel  is  less  uniform  and  varies  between 
0.50  and  0.79  A  cm-2.  The  current  density  for  the  TPB  under¬ 
neath  the  rib  (x>3mm)  also  decreases  rapidly  with  increasing  x 
and  is  determined  by  the  corresponding  TPB  oxygen  concentra¬ 
tion. 

It  is  worthy  noting  that  the  hydrogen  concentrations  at  the 
anode  TPB  are  quite  similar  in  both  designs.  The  effect  of  the  rib 
presence  in  the  anode  side  on  the  TPB  current  generation  is  com¬ 
pletely  dominated  by  the  effect  of  the  rib  in  the  cathode  side.  As 
far  as  the  hydrogen  transport  and  the  TPB  current  generation  are 
concerned,  the  presence  or  absence  of  the  anode  rib  is  inconse¬ 
quential.  As  the  average  electrical  path  for  the  cross-flow  design  is 
the  same  as  the  co-  or  counter-flow  design,  the  numerical  results  in 
the  present  work  may  be  also  applicable  to  the  cross-flow  design. 

3.2.  Optimizations  of  the  rib  widths 

Fig.  7a  shows  a  contour  map  of  the  average  cell  current  density 
as  a  function  of  the  rib  width  and  channel  width  for  the  reference 
cell.  Generally,  the  smaller  the  pitch  width  (dpitch  =  ^channel  +  drjb, 
Fig.  7a  is  only  for  dpitch<4mm),  the  higher  is  the  output  cur¬ 
rent  density.  For  a  given  pitch  width,  there  is  an  optimal  rib 
width  for  producing  the  maximum  cell  current  density,  as  can 
be  seen  in  Fig.  7a.  The  maxima  of  current  densities  increase 
with  the  decrease  of  the  pitch  width.  Flowever,  the  pitch  size 
may  be  limited  by  the  manufacturing  and  engineering  consid¬ 
erations.  For  a  given  pitch  width  as  indicated  by  the  dashed 
line  in  Fig.  7a  (dpitch  =  dchannei +  drib  =  constant),  the  output  cur¬ 
rent  density  also  depends  strongly  on  the  rib  width  and  a 
suitable  choice  of  the  rib  width  is  very  important  for  the  high 


Channel  0™)  ^ 


Fig.  7.  The  relationship  between  the  cell  output  current  density  and  the  channel  and  rib  widths:  (a)  the  output  current  density  vs.  the  channel  width  (dchan nel)  anc*  the  rib 
width  (dri5):  (b)  the  channel  widths  for  the  maximum  cell  outputs  vs.  the  pitch  widths. 
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Fig.  8.  The  dependence  of  the  optimal  rib  width  and  stack  cell  output  on  the  pitch 
width  for  ASRCOntact  varying  from  10-5  to  0.1  Q,  cm2 :  (a)  the  optimized  rib  width  vs. 
the  pitch  width  and  (b)  the  maximum  output  current  density  vs.  the  pitch  width. 


performance  of  a  SOFC  stack.  The  optimal  rib  width  was  found 
by  varying  the  rib  width  for  a  given  pitch  width  to  produce  the 
maximum  output  current  density,  which  is  the  maximum  on  the 
corresponding  dashed  line.  Fig.  7b  shows  the  relationship  between 
the  given  pitch  width  and  the  optimal  rib  width  that  produces 
the  maximum  cell  output.  It  is  interesting  to  note  that  the  opti¬ 
mal  rib  width  is  linear  to  the  pitch  width  and  can  be  expressed 
as  drib  =  0.075  mm  +  0.39dpitch  (0.5  mm  <  dpitch  <  5  mm)  for  the  ref¬ 
erence  cell  with  a  contact  resistance  of  0.05  £2  cm2. 

3.2.1.  The  effect  ofASRcontact  on  the  rib  optimization 

The  optimal  rib  width  is  determined  by  the  interplay  between 
the  gas  transport  and  the  current  conduction.  The  higher  the  rib- 
electrode  contact  resistance,  the  higher  the  rib  width  is  expected. 
Fig.  8a  shows  the  relationship  between  the  optimal  rib  width 
and  the  designated  pitch  width  for  a  range  of  ASRCOntact  varying 
from  10-5  to  0.1  £2  cm2.  To  be  more  representative,  the  fuel  with 
a  medium  hydrogen  molar  fraction  of  0.6  is  used  here.  As  can  be 
seen  in  Fig.  8a,  the  optimal  rib  width  shows  a  high  degree  of  lin¬ 
earity  with  the  pitch  width  for  any  of  the  given  ASRCOntact,  with  the 
slope  for  the  linearity  decreases  with  the  increase  of  ASRcontact-  The 
corresponding  current  densities  are  also  shown  in  Fig.  8b.  For  the 
small  ASRcontact  of  10-5  £2  cm2  and  the  small  pitch  width  of  0.5  mm, 


the  output  current  density  is  essentially  the  same  as  that  for  a  sin¬ 
gle  cell  (0.97  A  cm-2)  due  to  the  fact  that  the  optimal  rib  width  is 
very  small  and  the  rib  effects  on  ohmic  overpotential  and  concen¬ 
tration  overpotential  are  negligible.  However,  the  rib  effect  cannot 
be  ignored  even  for  such  small  ASRcontact  when  the  pitch  width  is 
large.  The  output  current  density  is  reduced  by  over  40%  for  the 
pitch  width  of  5.0  mm.  The  performance  degeneration  is  mainly 
due  to  the  ohm  loss  of  the  cathode  current  conduction  as  shown  in 
Fig.  5  and  discussed  above. 

As  shown  in  Fig.  8b,  the  output  current  density  of  the  stack  cell  is 
severely  influenced  by  the  ASRcontact-  Even  for  small  pitch  width  of 
1  mm  and  moderate  ASRcontact  of  0.025  and  0.05  £2  cm2,  the  optimal 
current  densities  are  only  56%  and  45%  of  that  for  the  single  cell, 
respectively.  It  is  no  surprise  that  the  best  reported  stack  cell  output 
is  less  than  half  of  that  for  the  best  single  cell  result. 

3.2.2.  The  effect  of  hydrogen  molar  fraction  on  the  rib 
optimization 

As  the  cell  may  be  fed  with  fuels  with  different  hydrogen  con¬ 
centrations  and  the  hydrogen  concentration  for  a  given  fuel  varies 
for  different  channel  locations  due  to  fuel  consumption,  the  sen¬ 
sitivity  test  of  the  optimized  rib  width  on  the  fuel  composition  is 
important  for  the  rib  design.  The  effects  of  inlet  hydrogen  molar 
fractions  on  the  optimized  rib  widths  and  the  output  current  den¬ 
sities  vs.  the  pitch  widths  are  shown  in  Figs.  9  and  10  for  ASRcontact 
of  0.05  £2  cm2.  Even  though  the  stack  cell  output  depends  on  the 
inlet  hydrogen  molar  fraction,  the  influence  of  the  inlet  hydrogen 
molar  fraction  on  the  optimal  rib  width  is  quite  limited.  The  results 
for  ASRcontact  of  0.025  and  0.1  £2  cm2  (not  shown)  provide  the  same 
conclusion.  The  optimal  rib  width  is  linear  to  the  pitch  width  for 
any  hydrogen  molar  fraction.  For  a  given  ASRcontact,  the  optimized 
rib  widths  are  essentially  the  same  for  any  hydrogen  concentration 
when  the  pitch  width  is  small.  Even  though  the  output  current  may 
vary  substantially,  the  ratio  of  the  optimal  rib  width  to  the  pitch 
width  changes  by  no  more  than  7%  for  any  hydrogen  concentration 
when  taking  the  hydrogen  molar  fraction  of  0.6  as  the  reference.  The 
“insensitivity”  of  the  optimal  rib  width  to  the  H2  concentration  in 
fuel  is  in  fact  due  to  the  “sensitivity”  of  the  contact  ohmic  resistance 
overpotential  to  the  rib  width.  The  optimal  rib  size  is  due  to  the  bal¬ 
ance  between  the  current  collection  and  the  gas  transport.  Smaller 
current  output  in  general  results  in  smaller  optimal  rib  width,  as  can 
be  seen  in  Figs.  9  and  10.  In  one  hand,  the  contact  ohmic  overpoten¬ 
tial  is  naturally  sensitive  to  the  rib  width  for  high  contact  resistance 
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Fig.  9.  Optimized  rib  width  vs.  the  pitch  width  for  different  hydrogen  concentrations 
in  fuel  and  ASRcontact  =  0.05  £2  cm2. 
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Fig.  10.  The  maximal  output  current  density  vs.  the  pitch  width  for  different  hydro¬ 
gen  concentrations  in  fuel  and  ASRCOntact  =  0.05  £2  cm2. 


rib.  In  the  other  hand,  the  optimal  rib  widths  are  relatively  small 
for  low  contact  resistance  ribs.  A  small  reduction  of  the  rib  width 
results  in  a  relatively  large  increase  of  the  ohmic  overpotential  with 
a  small  reduction  of  the  transport  overpotential.  Consequently,  the 
optimal  rib  widths  for  high  and  low  hydrogen  concentration  fuels 
are  not  very  different.  This  is  fortunate  as  the  optimized  rib  design 
is  of  general  applicability  for  different  fuel  operations.  As  the  H2 
molar  fraction  of  0.6  is  roughly  the  average  H2  concentration  in  the 
cell  channel  for  SOFC  stacks  operating  with  fuel  of  high  H2  molar 
fraction  (say,  0.97)  and  desirable  fuel  utilization  (say,  75%),  the  H2 
molar  fraction  of  0.6  is  used  in  the  following  discussions  in  order 
to  be  more  representative. 

3.2.3.  The  effect  of  cathode  porosity  on  the  rib  optimization 

Fig.  11  compares  the  optimized  rib  widths  and  output  current 
densities  for  cathodes  with  two  different  porosities,  keeping  all 
other  material  property  parameters  to  be  the  same  as  the  reference 
cell.  Fligher  porosity  is  beneficial  for  the  gas  transport,  resulting  in 
larger  output  current  and  allowing  for  the  use  of  wider  ribs.  Except 
for  the  rib  with  unrealistically  low  contact  resistance,  however,  the 
optimal  rib  widths  are  essentially  the  same  for  the  two  porosities, 
as  can  be  seen  in  Fig.  11a.  The  insensitivity  of  the  optimal  rib  width 
is  due  to  the  dominant  sensitivity  of  the  ohmic  overpotential  to  the 
rib  width. 

3.2.4.  The  effect  of  cathode  layer  thickness  on  the  rib  optimization 

The  rib  width  are  also  optimized  for  different  cathode  layer 

thickness  varying  from  50  to  400  p,m.  Except  for  cases  with  very 
small  pitch  widths,  the  optimal  rib  width  is  essentially  indepen¬ 
dent  of  the  cathode  layer  thickness.  As  shown  in  Fig.  12a,  however, 
there  may  are  appreciable  variations  in  the  output  current  densi¬ 
ties  in  cases  of  large  pitch  widths.  Reasonably  thick  cathode  layer 
is  beneficial  for  the  oxygen  transport  to  the  TPB  region  underneath 
the  rib  area  and  helps  to  reduce  the  overall  cell  concentration  over¬ 
potential,  as  can  be  seen  by  comparing  Figs.  6a  and  12b.  Moreover, 
reasonably  thick  cathode  layer  is  also  beneficial  for  the  reduction 
of  the  cathode  ohmic  overpotential  by  increasing  the  cross-section 
area  for  the  current  to  pass  through  to  reach  the  interconnect  rib,  as 
seen  by  comparing  Figs.  5a  and  12c.  The  nearly  constant  rib  width  is 
a  result  of  the  balance  among  the  improved  gas  transport  for  wider 
rib,  the  increased  current  density  for  wider  rib  and  the  reduced 
cathode  ohmic  overpotential  for  narrower  rib. 


(mm) 


Fig.  11.  Effects  of  cathode  porosity  on  the  optimal  rib  width  and  the  cell  output:  (a) 
the  optimal  rib  width  vs.  the  pitch  width;  (b)  the  maximum  current  density  vs.  the 
pitch  width. 

3.2.5.  The  effect  of  cathode  conductivity  on  the  rib  optimization 
Due  to  the  influence  of  various  factors,  e.g.,  porosity,  conduct¬ 
ing  particle  size  and  volume  fraction,  different  layer  structures,  the 
effective  cathode  conductivity  may  vary  for  different  designs.  The 
effective  cathode  conductivities  for  practical  cells  [9,15]  are  simi¬ 
lar  to  our  reference  value,  differing  by  less  than  a  factor  of  two.  To 
examine  the  effects  of  the  cathode  conductivity  on  the  rib  design, 
four  more  conductivities  corresponding  to  1/5,  1/2,  2  and  5  times 
of  our  reference  value  are  used  in  our  simulations,  while  all  other 
parameters  are  set  at  the  reference  values.  Fig.  13  shows  the  results 
of  the  optimal  rib  widths  and  output  current  densities  for  these 
cathode  conductivities.  Except  the  case  for  the  lowest  conductiv¬ 
ity  and  large  pitch  width,  the  effect  of  the  cathode  conductivity 
on  the  optimal  rib  width  is  small.  The  results  may  be  attributed 
to  the  dominance  of  the  rib  contact  ohmic  overpotential  over  the 
cathode  ohmic  overpotential  in  typical  cells.  For  the  cathode  with 
very  low  conductivity,  lowering  the  cathode  ohmic  overpotential 
by  shortening  the  current  path  is  desired.  Consequently,  the  opti¬ 
mal  rib  width  for  the  wide  pitch  design  is  relatively  large  for  the 
cathode  with  very  low  conductivity  in  comparison  with  the  high 
conductivity  cathode.  Nevertheless,  such  low  conductivity  is  not 
typical  for  practical  cells.  That  is,  the  optimal  rib  width  is  practically 
independent  of  the  cathode  conductivity. 
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Fig.  12.  The  output  current  density  for  different  cathode  layer  thicknesses  (Lc)  with  the  optimal  rib  width:  (a)  the  maximum  output  current  density  vs.  the  pitch  width;  (b) 
distributions  of  the  TPB  current  density  and  02  molar  fraction  for  Lc  =  200  p,m  and  dpitch  =  2  mm;  (c)  the  current  density  (A  cm-2 )  and  current  flux  vector  in  the  electrodes  for 
Lc  =  200  |xm  and  dpitch  =  2  mm. 


3.2.6.  Parameters  for  the  optimal  rib  designs 

As  discussed  above,  the  optimal  rib  width  for  maximum  cell 
output  depends  linearly  on  the  pitch  width,  i.e. 

^rib  =  A  -\~  B  X  dpitch  (14) 

The  intercept,  A,  and  the  slope,  6,  are  weakly  dependent  on  the 
hydrogen  molar  fraction  in  the  fuel  and  are  essentially  indepen¬ 
dent  on  the  cathode  porosity,  the  cathode  layer  thickness  and  the 
cathode  conductivity.  A  and  B ,  however,  are  dependent  on  the  rib- 
electrode  contact  resistance  that  varies  with  materials  as  well  as 


dpitch  <mm) 


processing  processes.  As  an  easy  to  use  engineering  design  guid¬ 
ance,  Table  3  lists  the  parameters  A  and  B  based  on  the  hydrogen 
molar  fraction  of  0.6  in  the  fuel  for  a  range  of  ASRcontact  covering 
the  cases  for  basically  no  contact  resistance  to  fairly  high  contact 
resistance.  For  practical  SOFC  stacks,  ASRCOntact  is  expected  to  be 
between  0.025  and  0.075  £2  cm2  and  the  parameters  for  ASRCOntact 
of  0.05  £2  cm2  may  be  viewed  as  generally  representative  for  opti¬ 
mal  designs.  The  parameters  for  ASRCOntact  larger  than  0.1  £2  cm2 
can  in  principle  be  obtained  by  extrapolation,  but  stacks  with  such 
ASRcontact  may  be  of  limited  practical  use. 


^.ch  (mm) 


Fig.  13.  The  optimal  rib  width  and  output  current  density  for  four  cathode  conductivities  varying  from  1/5  to  5  times  of  the  reference  conductivity  (crc):  (a)  the  optimized 
rib  width  vs.  the  pitch  width;  (b)  the  maximum  output  current  density  vs.  the  pitch  width. 
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Fig.  14.  The  output  current  density  as  functions  of  the  pitch  width  and  cathode  thickness  for  difference  cathode  conductivities:  (a)  crc/5;  (b)  crc;  (c)  5<rc;  (d)  the  maximum 
output  current  density  vs.  the  pitch  width,  where  the  solid  and  open  symbols  are  for  the  optimized  cells  and  the  cells  with  the  cathode  layer  thickness  of  50  p.m,  respectively. 


3.3.  Effect  of  the  cathode  layer  thickness  on  the  fuel  cell 
performance 

For  anode  supported  SOFCs,  thin  cathode  layer  with  thickness 
of  about  50  pan  is  typical  for  sake  of  easy  oxygen  transport  [2]. 
The  concept  is  valid  only  when  metallic  mesh  with  good  contact  is 
used  as  the  current  collector  that  corresponds  to  negligible  contact 
resistance  and  very  small  pitch  width.  For  practical  interconnect 
material  in  SOFC  stack,  the  contact  resistance  is  substantial  and  the 
rib  width  cannot  be  made  too  small.  As  discussed  above,  reason¬ 
ably  thick  cathode  layer  is  beneficial  for  oxygen  transport  to  the  TPB 
underneath  the  rib  as  well  as  for  broadening  the  cross-section  for 
current  passage.  As  shown  in  Fig.  12a,  the  benefits  may  result  in  sub¬ 
stantially  higher  current  output  for  practical  pitch  width.  Naturally, 
the  two  benefits  should  be  balanced  with  the  oxygen  transport  to 
the  TPB  region  underneath  the  channel  that  prefers  thinner  cathode 
layers.  Therefore,  the  optimization  of  the  cathode  layer  thickness  is 
of  importance  for  the  SOFC  stack  design. 

Fig.  14  shows  the  output  current  density  as  functions  of  the  pitch 
width  and  cathode  thickness  for  three  cathode  conductivities  cor¬ 
responding  to  1  / 5, 1  and  5  times  of  the  reference  value,  while  other 
parameters  were  the  same  as  the  reference  cell  and  the  rib  width 
optimized.  The  results  corresponding  to  the  cathode  conductivities 
of  1  \2  and  2  times  the  reference  value  have  also  been  obtained  and 
show  a  similar  trend.  The  benefit  of  the  cathode  layer  thickness 


optimization  is  shown  by  comparing  the  outputs  of  the  optimized 
cell  and  the  cell  with  a  50  p,m  cathode  layer,  keeping  the  other  con¬ 
ditions  equal.  As  shown  in  Fig.  14,  the  cathode  thickness  for  the 
maximal  current  output  is  usually  around  150  p,m  or  more.  The 
optimization  of  the  cathode  layer  thickness  benefits  the  large  pitch 
width  design  more  than  the  small  pitch  width  design.  The  optimal 
cathode  layer  thickness  is  larger  for  larger  pitch  width.  The  biggest 
factor  influencing  the  optimal  cathode  thickness  seems  to  be  the 
pitch  width,  and  not  the  cathode  conductivity.  Flowever,  the  cath¬ 
ode  layer  thickness  should  generally  be  less  than  300  pum  even  for 
large  pitch  width  as  there  is  basically  no  gain  in  the  output  current 
density  by  using  thicker  cathode  layer. 

The  current  results  of  the  cathode  layer  thickness  optimization 
are  only  primary.  More  systematic  examination  of  various  influenc¬ 
ing  factors  combined  with  the  detailed  architecture  of  the  cathode 
(e.g.,  the  function  layer  and  the  current  collection  layer  and  their 
microstructures)  will  be  required. 

4.  Summary 

A  2D  numerical  model  coupling  the  gas  transport  in  porous 
media  and  electrical  conduction  process  is  presented  and  detailed 
discussions  about  the  effects  of  the  interconnect  rib  on  the  per¬ 
formance  of  the  planar  SOFC  stack  cell  are  given.  The  interconnect 
rib  affects  the  stack  cell  performance  mainly  through  the  intrinsic 
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ohmic  polarization  due  to  the  rib  surface  contact  resistance,  the 
increased  cathode  concentration  polarization  due  to  oxygen  deple¬ 
tion  in  the  area  underneath  the  rib  and  the  increased  cathode  ohmic 
polarization  due  to  unfavorable  current  distribution.  The  intercon¬ 
nect  rib  may  drastically  affect  the  stack  cell  performance  and  the 
rib  design  optimization  is  of  high  engineering  importance. 

Systematic  optimization  of  the  rib  size  for  the  maximum  stack 
cell  output  and  sensitivity  analysis  on  the  optimization  results  have 
been  carried  out.  The  optimal  rib  width  is  found  to  be  dependent 
linearly  on  the  pitch  width.  The  smaller  the  pitch  width,  the  higher 
the  optimal  stack  output,  and  it  is  recommended  to  use  the  small¬ 
est  possible  pitch  width  allowed  by  the  engineering  considerations. 
The  parameters  in  the  linear  relationship  between  the  optimal  rib 
width  and  the  pitch  width  are  dependent  on  the  area  specific  con¬ 
tact  resistance  of  the  rib,  but  are  practically  independent  of  the 
electrode  porosity,  layer  thickness  and  conductivity.  The  parame¬ 
ters  are  only  weakly  dependent  on  the  hydrogen  molar  fraction  in 
the  fuel  and  the  results  for  medium  hydrogen  concentration  are  sat¬ 
isfactorily  representative.  The  simple  linear  relationship  between 
the  optimal  rib  width  and  the  pitch  width  and  the  parameters  given 
for  a  range  of  practical  contact  resistance  provide  an  easy  to  use 
guidance  for  the  broad  SOFC  engineering  community. 

Primary  results  on  the  effects  of  the  cathode  layer  thickness  on 
the  SOFC  stack  performance  are  also  presented.  The  optimization 
of  the  cathode  layer  thickness  is  particularly  beneficial  for  large 
pitch  width  and/or  low  conductivity  cathode  materials.  The  optimal 
cathode  layer  thickness  is  generally  in  the  range  of  100-400  p,m. 
The  larger  the  pitch  width,  the  thicker  the  cathode  layer  should  be. 
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